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Abstract 

Using a self-consistent delta-hole model the pion propagation in isospin asymmetric 
nuclear matter is studied. In neutron-rich matter, corresponding to heavy nuclei, 
a significant difference in positive and negative pion light-cone distributions is ob- 
tained leading to a nuclear enhancement of up antiquark distribution compared 
to the down antiquark one. This means that the sea-quark asymmetry in the free 
nucleon cannot be extracted model independently from an experiment on a neutron- 
rich nucleus. 

Key words: Isospin-asymmetric medium; Drell-Yan scattering; Antiquark flavour 
asymmetry. 



1 Introduction 



The meson-cloud model plays an important role in dealing with non-perturbative 
Quantum Chromodynamics effects in the nucleon. It has been used by several 
groups to interpret momentum distributions of sea quarks in the nucleon [1-5], 
measured in deep inelastic scattering. The standard approach is based upon 
the Sullivan process [6] , in which the only essential parameter is the cut-off in 
the pion-nucleon-nucleon vertex. 

Originally the emphasis was mainly on the description of the isoscalar u + d 
distributions. More recently, since the observed violation of the Gottfried sum 
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rule, showing an excess of d over u in proton, the interest focussed on the 
properties of u{x) — d{x) in the nucleon, whose x dependence was measured 
recently [7]. The asymmetry in antiquark distributions has been interpreted 
mostly in terms of the pion-cloud model (a review is given in ref. [8], see also 
[9]), and also in a soliton model in the large- A^^c limit [10]. 

Since pion properties are strongly affected by the nuclear medium, the pion 
cloud plays also an important role in modelling nuclear effects on deep inelastic 
processes. It was used in the past in connection with the EMC effect [11-13], 
leading to a few percent enhancement of the structure-functions ratio around 
X = 0.2. Similarly, most approaches predicted a nuclear enhancement of the 
pion field and hence the u{x) and d{x) distributions, leading to a noticeable 
increase in the predicted Drell-Yan (DY) cross-section ratio. On the other hand 
experimentally the DY scattering on nuclear targets [14] did not show much 
evidence for enhancement over nucleon targets. Several papers have dealt with 
this discrepancy, pointing out different mechanisms leading to reconciliation 
with the measurements [4,5,15,16]. 

For the purpose of presenting our results we find it convenient to separate them 
into two aspects. The calculation of experimentally measured proton-nucleus 
to proton-deuteron DY cross-section ratio is our first aim. The isospin asym- 
metric medium affects the various isospin states of the nucleon's pion cloud 
differently, leading to an excess of up antiquarks over down ones (even if the 
distributions in free proton are identical) for nuclei with more neutrons than 
protons. Since the charges of up and down quarks are different, the effect also 
shows up in the nuclear DY scattering cross section. In addition to presenting 
the DY ratio, one can also investigate possible nuclear effects on the difference 
u{x) — d{x) for a nucleon in the medium. It has been pointed out by Kumano 
that in a neutron-excess nucleus there could be medium effects contributing 
to this difference, coming from flavor-asymmetric parton recombination in the 
small Bjorken-x region [17]. The mechanism he considered contributes only at 
relatively small x values [x < 0.1) and the effect is also rather small (2-10% 
of the observed asymmetric sea in the nucleon). 

In the present calculation we extend the effective field theory approach to 
compute pion properties in isospin-symmetric nuclear medium [18] to the case 
of arbitrary proton and neutron densities. Then we use the pion-cloud model 
[16] to compute the medium-modified quark and antiquark distributions and 
the proton-nucleus DY cross section. We assume that the medium effects solely 
originate from a modification of the pion cloud. Q The nucleons are treated in 



^ It is known [19] that the pion and other mesons cloud cannot account completely 
for the antiquark distribution of the nucleon. For example, gluon splitting gives a 
sizeable contribution at small x and large Q^, but this contribution is approximately 
flavor symmetric [19] and thus should not modify our results significantly. 
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mean field approximation, while the delta-isobar and the pion are dressed self- 
consistently. The effects of short-range baryon repulsion are included through 
Migdal's g' parameters {g'j^pfi g^A^ 9'aa)- Sensitivity of the pion distribution to 
the delta-hole self-energy was already recognized in ref. [11], thus we carried 
out the computation with the self-consistently determined in-medium delta- 
isobar spectral function [18]. 

We compute the pion light-cone momentum distributions separately for the 
three charge states as a function of the asymmetry parameter P = {N — Z)/A. 
For /3 > the 7r~ distribution exceeds that of tt*^, which in turn is larger 
than the 7r+ distribution, as expected on the basis of particle-hole self-energy 
relationships. 



2 Drell-Yan cross section 



The Drell-Yan cross section for the process p -\- A ^ /I'^ii X is given by 
(suppressing the dependence) 

= —-^1^^-^ Z]4[^/(^i)^/(^2) + qf(xi)qf{x2)]dxidx2, (1) 



where the sum is over all flavors, and xi, X2 are the longitudinal momentum 
fractions carried by quark of the beam and target nucleons, respectively. By 
a suitable selection of kinematics the values of xi,X2 can be deduced from 
experiment [14]. 

If we consider the region Xi > 0.3 when the antiquarks in the projectile play 
a negligible role, the ratio of the proton-nucleus to proton-deuteron DY cross- 
sections takes on the form 



^ 2 da{pA) _ un/a{x2) + dN/A{x2) 
A da{pd) Up{x2) + dp{x2) 

I fix,) ^^'^^^'^'^ -dN/A{x2) ^2) 
^ Up{x2) + dp{x2) 

where f{xi) = {Aup{xi) — dp{xi)) / {4:Up{xi) + dp{xi)) is close to unity. Here 
Up and dp are antiquark distributions in the free proton, while Un/a and rf^r/yi 
are the antiquark distributions per nucleon in the nucleus, differing from the 
free-nucleon distributions by the medium modified pion-cloud contribution. 
Denoting the latter by t^w^/^ and Sd^^/A the DY ratio becomes 
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RAd = 1 + —7 — r^— 77 — 7 {^u^/a{x2) + Sd^/A{x2) 

+ pf{xi) dp{x2) - Up{x2) + {5u^/a{x2) - 5d^iA{x2)) 1(5 } . (3) 

We see from the above expression that apart from nuclear effects, leading to 
nonzero Sut^/a and 5d,r/A, for /3 7^ there is a nucleonic one [4], stemming 
from the nonzero value of the antiquark-distribution difference dp — Up in 
the free proton. This underlines the necessity to use parton distributions in 
accordance with latest dp — Up observations (for discussion on this point see 
section IV). For the relatively small asymmetries of interest for stable nuclei 
the ratio {5ut^/a — ^d.,^iA)/(i appearing in the above expression is practically 
independent of /3, leading to a linear dependence of RAd on it- 

For the change of antiquark distributions due to pion-cloud modification we 
use the convolution formula [16] 

where 5f'f'{ii), given by 

5r;{y) = ^ - r'^^iy)) + ^ [r'"\y) - ,(4) 



represents the change of the pion light-cone-momcntum distribution per nu- 
cleon in the medium. /'^"/"M i^f^" Ip/a-^ and /'^"Z" (/'""/'P) denote the distri- 
bution of tt" per neutron (proton) in medium and in free space, respectively. 
They are discussed in the next section. 



3 Pions in isospin asymmetric nuclear medium 



We consider a model consisting of pions, nucleons and delta-isobars in an 
infinite, spatially uniform system at zero temperature. The proton and neu- 
tron densities are given through their chemical potentials ^p and The 
equilibrium conditions for nucleons, delta-isobars and pions imply that the 
chemical potential for the neutral pion is zero, while those of charged pi- 
ons are: ji-j^- — — /i7r+ = Hn — l^p- The antiparticles of nucleons and isobars 
are neglected, but a relativistic kinematics is used. The nucleons are fur- 
ther treated in the mean-field approximation, with momentum-independent 
mass {M^p — Mp, M^n — Mn) and energy shifts (cp, Cn) modelling their bind- 
ing, i.e. En{p) = ^JmI~^ + c„ for neutrons and Ep{p) = ^M^^ + p^ + Cp 
for protons. The Schwinger- Dyson equations without vertex corrections are 
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then solved self-consistcntly for the delta-isobar and pion [18]. The pion self- 
energy consists of the particle-hole and delta-hole contributions, with both 
imaginary and real parts taken into account, assuring correct analytical prop- 
erties. Short-range baryon repulsion is taken into account through Migdal's 
9'nn J 9n At q'aa parameters. The sum- rules for spectral functions of pions and 
deltas are checked and found to be satisfied to 1-2%. 

The pion light-cone momentum distributions are then calculated, using the 
in-medium pion self-energy. Direct calculation of the diagrams corresponding 
to the pion emitted by the in-medium nucleon proceeds analogously to ref. 
[16], but separately for the three charge states of the pion and taking into 
account different neutron and proton densities (and thus different M^p, M^^ 
effective masses and Cp, Cn energy shifts). 

It is well known that f{y) can be expressed compactly in terms of an integral 
over the spin-isospin response function or imaginary part of the pion prop- 
agator. It is also possible to compute it directly from summing contributing 
diagrams [16], which represent emission of a dressed pion by a nucleon. The 
latter procedure we prefer for numerical reasons, since it needs computation 
of the pion self-energy in smaller energy- momentum region. The schematic 
expression for the pion light-cone distribution per nucleon has the form 

where represents the sum of relevant diagrams, p the nucleon density, 

ko — E{p) — E{p'),k = p — p' and p {p') is the momentum of incoming 
(outgoing) nucleon. 

The full expression generalizing the case for isospin-symmetric medium from 
[16] reads (for the special case of tt"*" on proton): 





2 /If ^''p yPpp-Ps oo 



PFp 



2n 



X [ de^F^r^Nik)^ \D^+{k^M\ (5) 



where ppp is the Fermi momentum of the proton, M the nucleon's mass, = 
{M^p + M^^l2^ Q-kNN the n^NN couphng and F^iVAr(A;) the form factor, while 
6 is the angle between p± and p'j^. is the full pion propagator, corrected 
for the presence of four-fermion couplings through Migdal's g' parameters, 
given as D^^ in ref. [16], p'^^,^ = {2z'M4Ml + p^J^' - {z'^ + 1)M^J'/^ 
if the argument of the square root is positive, otherwise is zero, and z' = 
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{—My +P3 + E{p) — Cn)/M^,. The expression for tt^ is obtained by swapping 
the indices n and p and for neutral pions the total distribution is a sum of two 
terms corresponding to emission by proton or neutron. 

The difference in distributions for the various charge states of pion basically 

comes from two factors. One is the Pauli blocking of the outgoing nucleon, 
which in neutron-rich medium restricts emission of vr^ (from a proton, creating 
a neutron in the final state) more than the emission of tt" (since a proton 
appears in this case in the final state) . The other effect is the dressing of the 
pion propagator, in which the particle-hole and delta-hole self-energies enter. 
Since the dominant contribution comes from the particle-hole contribution 
for > Z (neutron density larger than the proton one) the propagator 
is affected more than the one of the tt"'" (more details are given in the next 
section) . 

While the delta has been shown to play an important role in the asymmetry 
for free nucleons, the medium effects are negligible [16] and not included here, 
which is partially a consequence of the use of a soft pion-nucleon-delta form- 
factor, as obtained from a fit to pion- nucleon scattering [18]. Since the isobar's 
contribution for /3 7^ to Sutt — Sd.„ is of the opposite sign compared to the 
contribution of Eq. (5), as discussed in the next section, at small x inclusion of 
this term might result in a small decrease of the calculated isospin-asymmetry 
effect. 



4 Results and discussion 

For numerical calculation we used the proton to neutron ratio of tungsten, for 
which there are measurements of the DY cross section [14]. The asymmetry 
parameter in this case is /3 = 0.196 and for the Fermi momenta of protons 
and neutrons we chose ppp = 238 MeV and ppn = 272 MeV, corresponding 
to total nucleon density slightly below the saturation density. To take into 
account the different binding of protons and neutrons for the energy shifts we 
take Cp — 40 MeV, Cn — 42 MeV, with the common effective mass = 
M*p = 0.85 GeV, thus assuring the correct asymmetry energy of 28 MeV. 
For the pion-nucleon-nucleon vertex we use a dipole form-factor with cut-off 
A = 1 GeV. We checked that varying the cut-off in the range 0.9 — 1.1 GeV 
does not change appreciably the medium effect on the pion. 

In Fig. 1 we present the pion light-cone-momentum distributions for nucleons 

in medium and in free space (upper four curves), as well as the excess pion 
distributions (lower three curves). The 7r+ (vr^) distributions are per proton 
(neutron), while that of tt" is per average nucleon in the medium. The neu- 
tral pions see little difference between an isospin-symmetric and an isospin- 
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Fig. 1. Pion light-cone-momcntum distributions. Full line is for free nu- 
cleon, long dashed for vr^, short-dashed for vr^, dot-dashed for vr'^, with 
9'nn ~ 0-6)5'jvA ~ S'aa ~ O-"^- '^^^ three lower curves show the excess distribu- 
tions with respect to the free nucleon. 

asymmetric miclcar medium, as long as the total micleon density is the same 
(actually, the different mean-field shifts for proton and neutron may lead to a 
very small effect). However, for the case of neutron excess, the 7r~ distribution 
per neutron is larger than the 7r+ distribution per proton. This comes from 
a difference in particle-hole self energies, which is mainly responsible for the 
light-cone-momentum distributions. It is easy to understand the difference if 
we look at the imaginary part of the self energy as a function of pion's energy. 
If the energy is positive, the pion can excite a neutron from the Fermi sea 
to become a proton (above the proton Fermi sea) if its charge is positive. A 
negative pion can make from a proton in the Fermi sea a neutron above the 
neutron Fermi sea. Since there are more neutrons than protons in the Fermi 
sea, the absolute value of the imaginary part of the tt^ self energy will be 
larger than that of tt". Bearing in mind that in expression (5) the pion en- 
ergy is negative (for the dominating part of the integraton region) and using 
the relation that \m.'n.T^+{uj^k) = lTaIlT^-{—uj, k) we arrive at the inequality 
JIT /n/A^y'j ^ j"^^ ^^^^{y), in accordance with the numerical computation. We 
mention that for the delta-hole self-energies the relationship between the 7r~ 
and TT"^ is the opposite to that of particle-hole one, i.e. for positive energy that 
of TT" is larger (in absolute value) than the self-energy of tt"*"; a consequence of 
different isospin factors in the pion-nucleon-delta vertex. We see that for the 
studied asymmetry there is a significant difference in the distributions of three 
charge states. The distribution in isospin symmetric medium is very close to 
the 7r° distribution. 
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Fig. 2. The antiquark-distribution difference, {5ut^/a — 5^7^/^)7/5 per nucleon in 
the medium. Fuh hne is for g'j^j^ = 0.6, g'j^^ = g'^^ = 0.3, dashed hne for 
= 0.5,5^^ = fiAA_= 0-3 and dash-dot Hne for = 0.5, c/^^ = g'^^ = 0.4. 
The dotted hne shows dp — Up for comparison. 

The valence distribution of negative pions contains up antiquarks and thus 
in neutron-rich matter they outnumber the down antiquarks present in the 
positive pions, due to nuclear effects on the pion clouds. A comparison of 
this effect to the recently measured [7] dp{x) — Up{x) difference for proton is 
presented in Fig. 2. It shows the quantity {Sun/A — Sdn/A)/P compared to the 
mentioned dp{x) — Up{x), to which it is added in the expression (3) for the 
DY ratio Raci- The quantity {Sut^/a — Sd^^/Aj/P does not change appreciably 
with P up to its value of 0.2, but it is sensitive to the g' parameters. Existing 
experimental and theoretical information [15] suggests values of 0.55 - 0.7 
for g'^^ and 0.3 - 0.4 for g'j^^ and g'/^^. These values give results consistent 
with observed DY scattering for the isospin-symmetric calculation [16], and 
we use these values also in the present case. Taking momentum dependent g' 
parameters (as suggested in ref. [15]) could change the details of our results. 
Exploration of this and other effects we leave for a future publication. 

The free nucleon parton distributions are taken from ref. [20]. They fit the 
dp — Up and dp/up of ref. [7] quite well up to x = 0.15. However, at larger x 
the measured difference dp — Up and especially the ratio dp/up are poorly fitted 
[7]. To correct for this discrepancy, which would cause a significant increase 
in calculated RAd^ we impose a constraint dp{x) = Up{x) for x > 0.3 (by using 
the arithmetic mean value) and interpolate linearly between x = 0.15, when 
the unmodified distributions are used, and the point x = 0.3. In this way, the 
employed distributions fit the measurements of ref. [7] very well. From Fig. 2 
we see that the two terms, whose sum appears in the square bracket of Eq. 
(3), are of comparable magnitude, i.e. the nuclear effect of the isovector part 
of the DY ratio plays as important role as the nucleon antiquark asymmetry. 
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Fig. 3. DY cross-section ratio. Full (dash-dot) line is for asymmetric nuclear medium 
with f3 = 0.196 and g'j^j^ = 0.6 {g'j^fj^ = 0.5), dotted line for symmetric matter with 
b'nn ~ 0-6; in both cases g'j^^ = g'^^ = 0.3. The dashed line is for (3 = 0.196, but 
without nuclear pion effects. Experimental results from ref. [14] for W are shown as 
points with error bars. 

Since the squared charge of up antiquarks is four times that of down anti- 
quarks, the former enter in the DY cross section with correspondingly larger 
weight. This implies an enhancement of the DY cross-section ratio for neutron- 
rich medium compared to the isospin symmetric case. In Fig. 3 the ratio of 
cross-sections (per nucleon) is shown as a function of X2, where both in the 
numerator and denominator integration over Xi is performed for Xi > X2 + 0.2, 
corresponding to the experimental situation of ref. [14]. The full line (dash- 
dot line) corresponds to P = 0.196 and g'j^j^ = 0.6 {g'j^fj^ = 0.5), dotted line is 
for (3 = 0, while the dashed line is for asymmetric medium (/? = 0.196), but 
without medium effect on the pion cloud. Measurements from ref. [14] for W 
are shown as points with error bars. The errors are too large for any definite 
statement to be made on the isospin-asymmetry effect. 

We mention that the difference for larger X2 values between the (3 = case of 
Fig. 3 in the present work and Fig. 8 of ref. [16] is a consequence of the use 
of a more realistic dipole pion-nucleon-nucleon form-factor in the present cal- 
culation, and to a smaller extent due to a different set of parton distributions 
[20]. 

We remark that the small pion excess probability found in the present work 
would also lead to a pion contribution to the EMC ratio for a; ~ 0.1 which 
is smaller than in some other approaches [13]. However, this ratio is mostly 
sensitive to the nucleon self-energy in the medium and the role of pions is 
difficult to isolate. 

We conclude: i) a simultaneous experiment on a neutron-rich nucleus and an 
N ^ Z nucleus with a 5% accuracy should in principle be able to isolate the 
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see up-down asymmetry (term proportional to P in Eq. (3)); ii) about 50% 
of the u — d difference in a nucleus comes from nuclear effects. Therefore the 
assumption that nuclear effects are negligible as in ref. [4] cannot be justified. 
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